Dynamics of gas bubble growth 



a^ : in a supersaturated solution with 

o 



Sievert's solubility law 



in 

(N 

G. Yu. Gor, A. E. Kuchma 

o 

g . June 25, 2009 

> 

■ Department of Theoretical Physics, Research Institute of Physics, 

I Saint-Petersburg State University, 

^ ■ 1 Ulyanovskaya Street, Petrodvorets, St. Petersburg, 198504, Russia 

Abstract 

m ! 

\^ . This paper presents a theoretical description of diffusion growth of a gas bubble after its 

• nucleation in supersaturated liquid solution. We study systems where gas molecules com- 

o 

\ pletely dissociate in the solvent into two parts, thus making Sievert's solubility law valid. 

•'^ • We show that the difference between Henry's and Sievert's laws for chemical equilibrium 

X 

' conditions causes the difference in bubble growth dynamics. Assuming that diffusion flux 

is steady we obtain a differential equation on bubble radius. Bubble dynamics equation is 
solved analytically for the case of homogeneous nucleation of a bubble, which takes place 
at a significant pressure drop. We also obtain conditions of diffusion flux steadiness. The 
fulfillment of these conditions is studied for the case of nucleation of water vapor bubbles 
in magmatic melts. 
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1 Introduction 

Nucleation and growth of gas bubble in gas-liquid solutions are highly widespread phenomena 
both in nature and in technological processes. Gas bubble growth plays an important role in 
the production of microcellular materials and foams, as well as in glass production.^ Besides, 
the growth of water vapor bubbles dissolved in magmatic melt is a process which leads to 
volcanic eruptions. 

This paper presents a theoretical description of diffusion growth of a gas bubble in liquid 
solution resulting from a pressure drop (i. e. instant creation of supersaturation). It has to be 
noted that growth regularities of a solitary new phase embryo are crucial for the description 
of the evolution of the whole process of the phase transition in supersaturated solution. 

The description of gas bubble growth in supersaturated solution traditionally employs two 
rough approximations:'^"'' the flux of the dissolved gas toward the bubble is assumed to be 
steady; and the consideration is hmited only to bubbles with such a large radius that it 
is possible to neglect for Laplace pressure in comparison with the external pressure of the 
solution. 

In the present paper, while exploiting the steady approximation, considering the bubble 
from the very moment of its nucleation, we take into account time- dependent Laplace pressure 
in the bubble. This approach makes both gas density in the bubble and the equilibrium con- 
centration of the dissolved gas at the surface of the bubble time-dependent as well (hereinafter 
we will refer to number density of molecules using term "concentration"). 

Gas bubble growth in a solution taking into account Laplace pressure was considered as 
early as in 1950 in the classical paper by Epstein and Plesset,^ where the authors obtained 
the equation for the bubble radius as a function of time. In order to correlate the equilibrium 
concentration of the dissolved gas and the solution pressure, Ref.^ presupposed the fulfillment 
of Henry's law: i. e. the proportional dependence between these two values. Indeed, Henry's 
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law is fulfilled in many systems, e. g. in case of the solution of CO2 in magmatic melt;^ 
however, it is not universal: one cannot apply it to the case of the solution of H2O vapor. ^'^ In 
the latter case, which is crucial for practical reasons, Sievert's law is observed: the equilibrium 
concentration of the dissolved gas is proportional to the square root of the solution pressure. 

In the current paper we will obtain the equation for the bubble radius as a function of time 
similarly to the equation obtained in Ref.^ for the case of Henry's law. We will show that the 
difference between Henry's and Sievert's laws for chemical equilibrium conditions at the bubble 
surface causes the difference in bubble growth dynamics. We will obtain conditions when the 
steady approximation for the diffusion flux of gas molecules toward the bubble is applicable. 
We will show that, as the radius of the bubble increases, the steady condition becomes stricter; 
and, consequently, as a rule, the steady regime of bubble growth gradually gives way to the 
nonsteady one. 

The case of homogeneous nucleation of the bubble, which takes place at the significant 
pressure drop, will be supplied with the analytical solution for the obtained bubble dynamics 
equation. 

In section [7] of the current paper we will analyze whether the steady approximation is 
applicable to the case of gas bubbles in magma described by Navon^ and Chernov et al.^ for 
large radius of a bubble (neglecting Laplace forces). 

2 Equilibrium concentration of dissolved gas: Henry's 
law and Sievert's law 

Let us consider a gas solution in liquid which previously was in equilibrium state at temperature 
T and pressure Pq. The concentration of the dissolved gas in the solution under such conditions 
will be denoted as no. 
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Then we instantly relieve the external pressure to value 11 in such a way that solution 
becomes supersaturated. The temperature and the volume of the solution remain the same, 
thus value Uq still serves for the dissolved gas concentration. 

It is more convenient to express the state of the solution in terms of dimensionless variables: 
supersaturation ( and gas solubility s 

^^no-n^^ (1) 

where is the equilibrium concentration of dissolved gas, k is the Boltzmann's constant. 

Some time after the pressure drop a gas bubble nucleates and starts growing regularly. We 
assume that the bubble is in mechanical equilibrium with the solution, and its dynamics is 
governed only by diffusion process. We will not discuss the settlement of the mechanical equi- 
librium between the bubble and the solution, assuming that this settlement occurs much faster 
than the settlement of chemical equilibrium, as in Ref.^ The time of the bubble nucleation is 
considered as t = 0. The radius of the bubble will be denoted as R. FoUowing^"^ we treat 
the liquid solvent as non- volatile (which is true for systems considered in section [7!); therefore, 
we assume that the bubble consists of gas only, but not of the solvent vapor. The gas in the 
bubble is considered to be ideal. 

When the bubble is studied after sufficiently long time since its nucleation, its radius can 
comply with the strong inequality 

R > 2a/n, (3) 

where a is the surface tension of the pure solvent (it is true while the solution is considered 
diluted). Eq. ([2]) makes it possible to neglect the influence of Laplace forces on the bubble 
growth. Therefore, we can write the following equations for the pressure in the bubble Pr and 
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for the equilibrium solution concentration near the surface of the bubble ur. 

Pr = n, (4) 

riR = rioo- (5) 

The subscription oo denotes that the equilibrium concentration tt-oo is related with the equi- 
librium near the flat surface of phase separation {R — > oo). As follows from Eq. (jl]), the gas 
concentration in the bubble Ug is constant. Using the ideal gas law we have 

When Eqs. Q, dl]), dS]) are fulfilled, the description of bubble dynamics is trivial for 
the case of steady-state diffusion and can be described analytically even for the non-steady 
case.^°'^^ 

From the moment of bubble nucleation and till Eq. ([3]) becomes valid, Laplace forces 
influence bubble growth. Thus quantities Pr and ur become radius- dependent (and, therefore, 
time-dependent). So, for Pr we now have 

^« = n + ^. (7) 

In order to write the equation for ur we need to know the solubility law. For the sim- 
plest case of Henry's law the equilibrium concentration of dissolved gas is proportional to the 
corresponding pressure 

^ = ^. (8) 

The study of bubble dynamics with consideration of Laplace forces for the case of Henry's law 
was discussed in^' . However, Henry's law is valid only in systems where there no molecules 
dissociation occurs during gas dissolution. In another important case, when the gas molecules 
completely dissociate in the solvent into two parts, Henry's law is replaced with so-called 



Sievert's law^^'^'' , and, therefore, Eq. ([8]) is replaced with the following one: 




Sievert's law is fulfilled for the case of water vapor dissolved in a silicate melt.^^ The study of 
such solutions is important both for glass production^ and for volcanic systems^' ^ . 

The replacement of Eq. ([8]) with Eq. ([9]) means that the boundary condition for the gas 
diffusion problem will be different for Henry's and Sievert's laws. The change of boundary 
condition, as we will see further, leads to the change in bubble dynamics. 



3 Bubble dynamics equation 

After the nucleation of the bubble, when its growth can be considered regular (i. e. the 
bubble can not be dissolved by fluctuations), its growth is governed by the diffusion flux of gas 
molecules into it. In this paper we will study the case when the diffusion flux can be assumed 
as steady. The conditions when such approximation is valid will be discussed in section |6l 

Considering the equality of gas concentration at the bubble surface and the equilibrium 
concentration ur and the equality of gas concentration far from the bubble and the initial 
concentration no, we can derive a simple expression for the steady flux density jr, 

Jd = D . (10) 

Here D is the diffusion coefficient of gas molecules in pure solvent (we assume that the solution 
is diluted). 

Now let us write the expression for the number of gas molecules in the bubble. Exploiting 
the ideal gas law and using Eq. ([7]) we have 



3kT 



^ 2a 
H + — 
R 



(11) 
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Differentiating Eq. (fTTl) . using denotation ([6]), we obtain 



where 



dN , ^odR 
dt ^ dt 



(12) 



(13) 



is the characteristic size of the bubble. 



The material balance between the dissolved gas and the gas in the growing bubble gives us 



the following equahty 



(14) 



Substituting expression for diffusion flux density ( ITOl) and the rate of change of the number of 
molecules Eq. (1121) into material balance equation (fT4|) and then, exploiting Eqs. ([2]) and ([6]), 



we have 



RR 



Ds 



no - riR 



(15) 



Using Eqs. ([I]), ([7]), ([9]) and (fT3|) . the fraction in the r. h. s. of Eq. ( IT5|) can be expressed as 



no - ur 



(16) 



Eq. (fT6|) allows us to rewrite the equation of bubble dynamics (fT5|) in its flnal form: 



RR 



R^ 
R 



Ds 



(C + 1) 



3^ 
2 R 



[17) 



It should be noted that the well-known paper by Cable and Frade presented only a special 
case of Eq. f|T7|) . when = — 1 (Eq. (34) in Ref.^) ; such value of supersaturation corresponds 
to the dissolution of gas bubble in the pure solvent. The form of this equation for the general 



case presented in the current paper has not been obtained previously. 
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4 Fluctuational nucleation, critical bubble and initial 
condition for the bubble growth 

All aforesaid (sections [H [21 [3] and the obtained equation for the bubble growth dynamics f|T7j) ) 
can be applied to the cases of both homogeneous and heterogeneous nucleation. Below we will 
consider only the homogeneous nucleation case. 

Homogeneous nucleation of a gas bubble in a supersaturation solution means fluctuational 
mechanism of its appearance and requires a significant pressure drop (of the order of 10'^ 
times) . It is these conditions of bubble nucleation that exist in magmatic melts during volcanic 
eruptions^ . Since we decided to describe only regular growth of a bubble, we need to exclude 
the very process of nucleation from our examination and consider a bubble when it is already 
supercritical. 

Under the notion of critical bubble we traditionally^' understand such a bubble which 
radius Rc corresponds to the extremum of work of bubble formation. A critical bubble is in 
mechanical equilibrium with solution at the initial pressure Pq in chemical equilibrium 
with solution with concentration hq. These two conditions together with the solubility law 
unambiguously define the value of Rc. 

The mechanical equilibrium condition, with Eq. ([7]), evidently gives us 



and chemical equilibrium condition allows us to relate value Pq to the supersaturation (. 
From Eq. iQ we have 




(19) 



or, using Eq. ([T]), 



Po = n(c + i)'. 



(20) 
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Finally, substituting Eq. fl20|) in Eq. f[T8|l for the radius of a critical bubble, we have 

= n[(c + i)2-i]- ^^^^ 

Using Eq. fl2T]) and Eq. f|T3|) we can write the relation between R^j and -Rc 

Ra = l [iC + ly - 1] Rc (22) 

which will be exploited later. 

Here it is important to make a remark considering quantity s - gas solubility. Eqs. and 
((2]) give us 

UnkT 1 , , 

If Henry's law was fulfilled, we would have C + 1 = Pq/U and, therefore, 

> = (24) 

s would be a tabular value defined only by the initial state of the solution (before the pressure 
drop). 



Since here we consider Sievert's law, we have (see Eq. ([20])) C + 1 = v^^o/n; thus instead 
of Eq. ([21D, we obtain 

s = (25) 

Here it is convenient to introduce constant K as the coefficient of proportionality in Sievert's 
law: K = Uq/ \/Po. It allows us to rewrite Eq. ( 125|) in the following form: 

Eq. (l26l) shows us that in the case of Sievert's law the solubility value s depends on the final 
state of the solution. In the current paper a solubility value is understood as gas solubility at 
the final pressure H - after the pressure drop. 



9 



To deal with the dynamics equation f|T7j) one needs to provide it with a reasonable initial 
condition, i. e. to choose the value of constant initial radius of the bubble, in the following 
equality 



In papers^' ^ there was no special meaning assigned to the value of Rf. the reason of the 
appearance of the bubble was totally excluded from discussion. Here we assume that the 
bubble nucleated fiuctuationally, i. e. crossed the barrier corresponding to radius Rc- It 
means that Ri has to exceed Rc, but evidently we cannot use the radius of a critical bubble 
Rc as the initial value for the radius. 

A bubble that nucleates fiuctuationally in the solution is capable of regular growth if it 
has passed and, moreover, moved away from the near-critical region where fluctuations are 
still strong enough. Thus we choose, following Ref.,^^ Ri = 2Rc. Such value guarantees the 
absence of fluctuations and, as we will see further, provides us with convenient expressions. 

The necessary condition for the fluctuational nucleation of a bubble is a significant pressure 
drop (-Po/n ~ 10'^) and, consequently, due to Eq. fl20|) . high supersaturation: C ~ 30 — 40. 
Further we will use the following strong inequality for supersaturation: 



5 Time dependence of bubble radius 

Let us now solve the differential equation f|T7|) for the time dependence of the bubble radius 
with the initial condition ^FJ\ of homogeneous nucleation of the bubble. 

The solution of this equation in general case is too cumbersome and even turns out to be 
unnecessary when inequality (!28|) is fulfilled, as we will ensure in the current section. Here we 



R{t)\t=o — Ri- 



(27) 



C»i. 



(28) 
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present solutions of this equation for two particular cases: 



(29) 



and 



R:^ Rc 



(30) 



When strong inequality ( !28l) is fulfilled, these two cases cover the whole range R > 2Rc of 
regular growth of the bubble radius, and that is the reason why general solution of Eq. (fTTI) 
is not necessary for the system under consideration. 

Let us begin with the case when inequality ( l29l) is fulfilled. In this case we can omit 1 in 
comparison with R^/R in the 1. h. s. of Eq. (fT7|) : and we can also omit 1 in comparison with 
fraction 3R^/2R in the r. h. s. of this equation. Thus we have 

Ds 



R 



Ra 



(C + 1) 



3 i?cr 

2~R 



(31) 



Separating variables and exploiting Eq. (l22l) for i?^- under the square root, we can rewrite Eq. 
(1311) in the form which make its possible to integrate it 

dR Ds{C + 1) 



, „ -dt. 

I-^RJR Ra 

Integrating Eq. (|32l) with initial condition ( !27l) we finally have 



(32) 



R-2Rc + RM 



R_ 

Rr 



1 \ +2jR, 



(yi - +Rc In 



V2 + lVR-^/R'c 
y/R+y/lTc 



Ra 



(C+iK 

(33) 



This formula is different from Eq. (3.4) in Ref.^^ , which means that, when inequality (l29l) is 
fulfilled, there is a significant difference in the character of growth between Sievert's law and 
Henry's law. 

Now let us proceed to the other case, Eq. fl30|) . At first, using Eq. fl22l) we can rewrite Eq. 
f[T7l) equivalently in the form of 

R[R + Ra] 



1 I Rc 



Ds{C+l). 



(34) 



(C+l)2 ^ R 
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This form makes it obvious that, when both strong inequahty fl3Ql) and inequahty fl28|) are 
fulfilled, the whole square root in the denominator of the 1. h. s. of Eq. fl34|) can be omitted 
in comparison with 1. Therefore, we have 

R[R + R^]=Ds{C + l). (35) 

This expression can be easily integrated. But initial condition (l27l) is not justified when 
strong inequality ( l30l) is fulfilled. There is arbitrariness in the choice of the initial condition 
for integration of Eq. (135!) . but the most convenient is to choose a condition at such a radius 
which simultaneously satisfies Eqs. and (pi . e. g. 

Rm^t^=R„, = V^Rjf.. (36) 

Due to Rm <^ -Rct; we can use Eq. (133|) to obtain the explicit value for time tm defined in Eq. 
(l36l) . Integrating Eq. ( 1351) with initial condition (l36l) . and then, substituting the expression 
for tm, we have 
R^ 

— — h RcrR + 2i?^^^i?j?^^ + RcrRc 

(37) 

Time dependences of bubble radius, when R <^ R^ (Eqs. ( l33l) and ( 1371) ). are presented 
graphically in Fig. [T]for Pq/H = 10"^. Bubble radius R is measured in units of Rc and time t 
is measured in units of tttt^- 

With the increase of R the contribution of the third and fourth addends in the 1. h. s. of 
Eq. (137|) decreases; and ai R = R^ we can already write 

— + R^R = DsCt. (38) 

In ( 138|1 we also took into account strong inequality (|28|) . Eq. (|38|1 is identical to Eq. (3.11) in 
Ref.^^ . It means that, when inequality R > R„ is fulfilled, any difference in the character of 
growth between Sievert's law and Henry's law disappears. 
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5 
2 



23/2 



Ds{C + l)t. 



When the bubble radius becomes as large as 

R:^ Ra, (39) 
Eq. (l38l) transforms into the well-known Scriven's^'' self-similar dependence for the steady-state 



case 



R^ = 2DsCt. (40) 

This trend of Eq. fl38|) toward Eq. fj40l) was discussed in detail in Ref.^^ This trend is also 
presented graphically in Fig. [2] for Pq/H = 10'^. Bubble radius R is measured in units of Rc 
and time t is measured in units of 



Ds(C+l)- 



6 Conditions for the diffusion flux steadiness 

The diffusion flux of molecules toward a growing bubble can be considered steady when the 
bubble growth is slow enough in comparison with the "diffusion cloud" growth. To be more 
exact, the radius of the bubble has to be much smaller than the radius of this cloud - the 
diffusion length. We can express it as 

R<^{DtRy/\ (41) 

where is the characteristic time of the bubble radius change, tji = R/R is the time in which 
the bubble radius changes significantly. Evidently, Eq. (14T|) can be rewritten as 

RR/D]'^ <^l. (42) 



We can make this condition more explicit by means of Eq. ([T 

It should be noted that condition ( l43l) is valid regardless of whether the bubble nucleation 
was homogeneous or not (the same can be applied to Eq. ( |T7I) ). Graphical dependence of 
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parameter characterizing the steadiness of bubble growth (L h. s. of Eq. fH3|l ) on bubble 
radius R is presented in Fig. [31 In Fig. [3] bubble radius is measured in units of Rc and 
presented in logarithmic scale, Pq/TI = 10^, solid curve: s = 1, dotted curve: s = 10~^. 

Now, exploiting Eq. (143!) . let us obtain conditions of steadiness of the diffusion flux at 
different characteristic sizes that we have for the case of homogeneous nucleation of the bubble. 
Obviously, the larger the bubble is, the stricter the condition for steadiness becomes (see Fig. 
[3]). The general condition of the steadiness of bubble growth at any time is the condition at 
i? > i?^. From Eq. (jl3D we have 

/i\ 1/2 

s'/' « ( ^ . (44) 



As we will find in the following section, Eq. (jH]) is difficult to satisfy. But it is evident that in 
order to satisfy this condition one needs to decrease the solution supersaturation ( (the case 
of heterogeneous nucleation) or, leaving supersaturation ( constant, decrease the gas solubility 
s. From Eq. (120|) for the supersaturation we have 

C=y§-1. (45) 
Then, using Eqs. fHSjl . let us rewrite product s( in the following form 

Eq. P6|) shows us that in order to weaken limitation fjUl) one needs to increase final pressure 
n, leaving ratio Pq/H constant. 

The condition of steadiness at the very beginning of the regular growth of the bubble, can 
be found, exploiting Eq. (143|) with R = 2Rc. 

s'/^ « C'/'. (47) 

1 /2 

While deriving Eq. (jlT]) . the numerical coefficient (3 (l — l/\/2)) in its 1. h. s. was replaced 
with 1 for shortness. 
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7 Bubble growth in magma during volcano eruptions 

This section contains the analysis of the question whether the steady growth conditions ob- 
tained in section [6] are fulfilled for bubble nucleation in volcanic systems. Previously^' ^ the 
study of such systems exploited steady approximation without the analysis of its applicability. 

It should be noted that usually there are several different gases dissolved in magma (e. g. 
HCl, HF, H2S, SO2, H2 and CO2), but the principal component is still water vapor^^ . For 
the analysis of steady conditions discussed below the presence of other gases is not important; 
therefore, we will consider the case of homogeneous nucleation of water vapor bubbles in 
magma. 

It is well-known that Sievert's law is fulfilled for gases where diatomic molecules dissociate 
during gas dissolution^^ . While water has a triatomic molecule, its vapor in magma obeys 
Sievert's law. This happens because each water molecule together with oxygen from silicate 
melt forms two hydroxyl groups (see Eq. (2) in Ref.^^): 

H2O + = 20H. (48) 

For our evaluations we will use parameters from Ref.^ ; we have 

Po = 100 MPa n = 0, 1 MPa T = 1150 K w = 3% p„, = 2300 kg/m^ (49) 

Here w is gas mass fraction of the dissolved gas (water vapor), and Pm is magma density. 
Notwithstanding the fact that liquid magmatic melt is a solution of several gases, none of 
these gases are surfactants and the solution is diluted; therefore, we can characterize such 
solution with surface tension a. Navon et al.^ gives the following values for the surface tension 
of magmatic melts: a ~ (0.05 — 0.1) N/m. 

Let us express the values of s and ( using data given. From Eq. fHSl) . in accordance with 
Eq. f l49p . we have ( ~ 31. Then we need to calculate Uq and substitute it into Eq. 0251] 
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to obtain the value of solubility. As long as we are given mass density of magma and mass 
fraction of gas, it is convenient to write 

'^o = Po , (50) 

P 

where po is mass density of the dissolved gas, Na = 6 x 10^^ mol~^ is the Avogadro constant 
and /i = 1.8 X 10"^ kg/mol is the molar mass of the dissolved gas (water). Finally, we need to 
express mass density of the dissolved gas. Evidently, we have 

Po = w^Pm, (51) 

and, therefore, 

NAkT 

Using data (H9ll in Eq. (l52l) . we have s ~ 12. 

Now we can see that for volcanic systems, where the pressure drop is of the order of 10'^ and 
solubility is more than 1, condition (jH]) is violated; and steady approximation is not valid for 
radii of the order of Ra- Even in the very beginning of bubble regular growth, when R = 2Rc, 
the steady condition (1471) is fulfilled only at its breaking point: the values of C^^^ exceed the 
value of s^/^, but these values are of the same order of magnitude. 

Let us make the following interesting observation. Since for R^ R^r the bubble dynamics 
is identical for both Henry's and Sievert's solubility laws, the condition of steady growth for 
-R 3> -Rcr is also the same (see Eq. flH|) above and Eq. (2.14) in Ref.^^). For the case 
of homogeneous nucleation, when the pressure drop Pq/H ~ 10'^, the steady condition at 
R:^ Ra rule, violated in cases of both Henry's and Sievert's laws. For Henry's law, 

when s ~ 10~^, it is violated due to high supersaturation values ( ~ 10^. For Sievert's law, 
when, due to Pq/H ~ 10^ and Eq. fl20|) . corresponding supersaturation values are significantly 
less C ~ 30 — 40, values of gas solubility are significantly higher than for Henry's law, and that 
is the reason of the violation of the steady condition. 
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8 Conclusions 

In the presented paper we obtained the equation to a bubble growth dynamics in the super- 
saturated gas solution with Sievert's solubility law. We showed that, as long as the bubble 
size is such that Laplace forces influence the bubble growth, the difference between Henry's 
and Sievert's laws in terms of chemical equilibrium conditions results in the difference in bub- 
ble growth dynamics. When Laplace forces become weaker, the bubble growth dynamics is 
identical in the cases of both Henry's and Sievert's laws. 

While obtaining the dynamics equation we assumed the diffusion fiux to be steady. We 
obtained conditions when this steady approximation is applicable. We showed that usually, as 
the radius of the bubble increases, the steady regime of bubble growth gradually gives way to 
the nonsteady one. Application of the obtained conditions for the volcanic system consisting 
of water vapor dissolved in silicate melt showed that the process in such system cannot be 
considered as steady. 
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Figure 1: Time dependence of bubble radius R from Eq. ( j33l) (dotted curve) and Eq. (j37j) 
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Figure 2: Time dependence of bubble radius R from Eq. fj40l) (dotted curve) and Eq. ([3 

(solid curve). Time t is measured in units of jj^^- 
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